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Abstract. We propose a model to describe a gas of pyramidal molecules 

interacting via dipole-dipole interactions. A cooperative effect induced 

^ , by the interaction modifies the tunneling properties between the classical 

(~| ' equilibrium configurations of the single molecule. The model suggests 

that, for sufficiently high gas density, the molecules become localized in 
these classical configurations. On this basis it is possible to explain the 
shift and the disappearance of the inversion line observed upon increase 

fj \ of the pressure in a gas of ammonia or deuterated ammonia. The same 

mechanism also accounts for the presence of stable optical activity of 
certain pyramidal molecules. We discuss the concept of environment 
[v> , induced superselection rule which has been invoked in connection with 

S_j ■ this problem. 



1 Introduction 

In this work we consider the behavior of gases of pyramidal molecules, i.e. 
molecules of the kind XY^ like for example ammonia NH^, deuterated ammonia 
ND^, phosphine PH^ and arsine AsH^- We begin by presenting an outline of the 
physical problem and of its history that dates back to the early developments of 
quantum mechanics. 
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For the single pyramidal molecule, owing to the great differences of charac- 
teristic energies and times, we can perform some adiabatic approximations. The 
electronic motion can be separated from the nuclear one (Born-Oppenheimer ap- 
proximation) . We can further separate the one-dimensional motion of inversion of 
the nucleus X across the plane containing the three nuclei Y . The form of the 
effective potential for this motion is a double well which is symmetric with respect 
to the inversion plane R, . The corresponding non degenerate eigenstates must be 
symmetric or anti-symmetric with respect to the inversion plane. Due to tunneling 
across the finite potential barrier, these eigenstates are delocalized in the two min- 
ima and, for energies below the barrier height, are grouped in doublets, i.e. couples 
of states with a relative splitting in energy small in comparison with the distance 
from the rest of the spectrum. For the pyramidal molecules under consideration, 
the thermal energy fc^T at room temperature is much smaller than the distance 
between the first and the second doublet so that we can reduce the problem to the 
study of a two-level system corresponding to the symmetric and anti-symmetric 
states of the first doublet. 

The existence of delocalized stationary states is clearly in disagreement with the 
usual chemical view which, relying upon classical theory, considers the molecules 
as objects with a well defined spatial structure. In particular, for the molecules 
under consideration the classical view predicts one of the two pyramidal configura- 
tions corresponding to the nucleus X localized in one of the wells of the inversion 
potential. 

The quantum prediction of stationary delocalized states implies the presence 
of a line in the absorption spectrum, the so called inversion line, at a frequency 
V = h~^AE, where AE is the energy splitting of the first doublet. The experi- 
mental results are strongly dependent on the different kind of atoms that form the 
pyramidal molecule. Experiments performed with NH3 and ND^ reveal the exis- 
tence of an inversion line, while this line is not found in the case of AsH^ and PH^. 
These results cannot be considered a demonstration of the existence of degenerate 
classical states for arsine and phosphinc. The values of the energy splitting for these 
molecules is so small, see Table W, that the inversion line, even if it existed, would 
not be detectable. However, there is a stronger evidence suggesting the existence 
of stationary localized states, at least for arsine (similar experiments are, as far 
as we know, not available in the case of phosphine). It has been found that some 
molecules of the kind AsYWZ, i.e. molecules obtained by substituting hydrogens 
with different atoms, show stable optical activity |g|. This means that the station- 
ary states of the AsHWZ molecule are chiral, i.e. they cannot be superimposed on 
their mirror image by a rotation, in agreement with the classical prediction of lo- 
calized pyramidal states. The discrepancy between the presence of optical activity 
for some kind of molecules and the quantum prediction of delocalized stationary 
states was pointed out for the first time in |4| and is usually referred to as Hund 
paradox or chirality problem. 

The frequency P of the inversion line measured in NH3 and ND3 gases shows 
a dependence on the gas pressure P ||, |[ 0. Starting from the expected value 
h~^AE at P ~ 0, PiP) decreases by increasing P and vanishes at a critical pressure 
Per — 2 atm for NH^ and Per ~ 0.1 atm for ND^. No quantitative theory has 
been proposed so far for this phenomenon. 

As early as 1949, in a short qualitative paper H Anderson made the hypothesis 
that dipole-dipole interaction may induce a localization of the molecular states. In 
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NH^ 


ND^ 


PHs 


AsHs 


AE (cm^^) 


0.81 


0.53 X IQ-^ 


3.34 X 10-^4 


2.65 X 10~^« 


^ (Debye) 


1.47 


1.47 


0.58 


0.20 



Table 1 Energy splitting of the first doublet AE and electric dipole moment 
fj, for the molecules NH3, ND3, PHs,, AsHs B. We remind that 1 cm"'^ = 
1.984 X 10-23 J and 1 Debye = 3.33 X lO-^*" Cm. 



this way the important idea was introduced that inter-molecular interactions may 
be responsible for the observed phenomena. Thirty years later Pfeifer studied a 
two-level system modeling a molecule which has an almost degenerate ground state 
and is coupled to the quantized electromagnetic field ^. He argued that the cou- 
pling of the molecule to the radiation field yields two symmetry-breaking effective 
ground states belonging to two different sectors of the Hilbert space separated by 
a superselection rule. These states are localized and one is the mirror image of 
the other. Later it was argued that this localization cannot take place at finite 
temperature [|o[ |l], |l|. 

A quantitative discussion of the effects induced by coupling a single molecule to 
the environment constituted by the other molecules of the gas was made in [|3[ |lj| . 
In this work it was shown that, due to the instability of tunneling under weak 
perturbations ||l^, |l^ , the order of magnitude of the molecular dipole-dipole inter- 
action may account for localized ground states in the case of arsine and phosphine 
while delocalized states were predicted for ammonia. Therefore, the appearance 
of chirality was interpreted as a collective effect in a system of infinitely many de- 
grees of freedom. This suggests that a phase transition may be invoked to explain 
the behavior of different pyramidal molecules or the behavior of the same kind of 
molecule under variation of thermodynamic parameters like pressure. 

We have implemented this idea by constructing a simplified model of a gas 
of pyramidal molecules which exhibits the desired properties and allows a direct 
comparison with experimental data. Our model predicts, for sufficiently high inter- 
molecular interactions, the presence of two degenerate ground states corresponding 
to the different localizations of the molecules. This means that there is a quantum 
phase transition for a critical value of density (or pressure in the case of constant 
temperature). A study at non zero temperature shows that this transition affects 
the behaviour of the gas at room temperature and gives a reasonable explanation of 
the experimental results ^, ^, 0|. In particular, it provides a correct description of 
the disappearance of the inversion line of NH^ and ND^ on increasing the pressure. 

In the present paper, after a description of the model, we discuss its mean field 
approximation and the appearance of a degeneracy for the ground state. We then 
analyze the linear response to an external field at zero temperature to obtain the 
low energy spectrum and its dependence on the inter-molecular interaction. This 
shows that the frequency of the inversion line decreases upon increasing density 
and vanishes at the critical density for which the ground state becomes degenerate. 
For a treatment at finite temperature and a detailed comparison with experiments 
we refer to our paper pT) ]. 

In the last section we discuss the interpretation of the chirality problem in 
terms of a superselection rule. 
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2 Description of the model 

As we mentioned in tlic introduction, for an isolated molecule XY^ we can 
adiabatically separate the electronic motions from the nuclear ones. The latter es- 
sentially reduce to the rotational motion of the molecule, the vibrations of the nuclei 
around their equilibrium positions and the inversion motion of the nucleus X across 
the plane of the nuclei Y . At room temperature the rotational and vibrational de- 
grees of freedom are much faster then the inversion motion and can be taken into 
account through a renormalization of the potential barrier for the inversion motion. 
We can further reduce the inversion dynamics to that of a two-level system corre- 
sponding to the ground and first excited eigenstates of the corresponding symmetric 
double well potential. For a general reference see (l) . 

Consider now a gas of interacting molecules. If the motion of the centers of 
mass of the molecules is slow in comparison with the two-level dynamics of each 
isolated molecule, the translational degrees of freedom can be considered frozen. In 
this case we can reduce our model to a lattice of unit length (. = q~^/^, where g is the 
gas density, with a two-state system at each site. By evaluating the inversion time 
hAE~^ with the data of Table |l|, we see that the lattice approximation is reasonable 
in the case of NH3 and ND3 up to room temperature. This approximation has been 
discussed in general ||l8| and we refer to this paper also for references to previous 
related works. 

In order to define a suitable Hamiltonian for our lattice model, let us begin 
with the isolated two-level system. We choose the Hamiltonian operator on the 
two-state Hilbert space to be 



2 '^' 






(2A) 


where a^ is the Pauli matrix 














(2.2) 


The corresponding eigenvectors are 








">^;^(0 






(2.3) 


'^>^;^(-0' 






(2.4) 


with eigenvalues —AE/2 and +AE/2, respectively, 
systems, the Hamiltonian is 


For 


N independent 


two-state 


^ AE 
i?o = > , ^li (8) I2 «)...«' (T- 


'(E,. 


■ • Iat, 


(2.5) 



where the subscript indicates the single particle space to which the operators are 
applied. 

In order to modeling the interaction among the molecules we briefly recall the 
discussion made in [Q. If we consider the whole gas it is possible that, owing to 
some fluctuation in the state of the system, a single molecule becomes localized. 
This molecule consequently acquires an electric dipole moment (i that polarizes the 
surrounding gas. The latter, in its turn, produces an electric reaction field Ef/, 
collinear with fi, so that the interaction fi ■ E^ tends to stabilize the localized 
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state of the initial molecule. In other words, there appears a cooperative effect 
which tends to stabilize the localization. To mimic this effect in the Hilbert space 
where Hq acts, we choose the interaction Hamiltonian 

N N 

Hint ^^ Yl 3{hj) li «)...«) erf (8)...(8) cr| «)...«) Itv, (2.6) 

where g{i,j) — g{j, i) < and a^ is the Pauli matrix 

whose eigenvectors are the localized states 

\L) = (^ J)=-^(|1) + |2)) (2.8) 

\R) - ( 5)=^(|1>-|2)), (2.9) 

with eigenvalues +1 and —1, respectively. The total Hamiltonian 

H = Ho+H,nt (2.10) 

contains two competing terms with Hi„t which tends to favor stationary localized 
states while Hq favors delocalization. According to the classical dipole-dipole in- 
teraction formula, we expect g{i,j) to be proportional to the square of the electric 
dipolc moment /_i of the molecules and to decrease with the cube of the distance 

l.9(M-)l « ^J^. (2.11) 

In the context of a mean field approximation, it is useful to introduce the 
parameter 

G^-J29ihJ), (2.12) 



whcr Nc is a cut-off. By neglecting border effects and using (2.11), we write 

G^C^, (2.13) 

where C is a positive constant which will be considered as a phenomenological 
parameter to be deduced, in the case of NH^ and ND^, from comparison with the 
experimental data on the inversion line. 

3 Ground state in the mean field approximation 

The mean field approximation for the ground state corresponds to minimizing 
the energy functional {ip \ H \ ip) within the class of states of the form 

|V') = |A)...|A), (3.1) 

where 

|A) =a|l)+6|2) (3.2) 

and 

\a\' + \b\' = l. (3.3) 
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This problem reduces to the solution of the non linear eigenvalue equation 

MA)|A)=r;(A)|A), (3.4) 

where h{X) is the self-consistent single particle Hamiltonian 

A F 
h{X) = a^- -Ga^\\a'\\). (3.5) 



The non linear term in (3^) accounts for the interaction of a single molecule with 



the reaction field and represents the polarization energy of the gas. 



We find |A) by solving Eq. (3^) and evaluate the energy functional {ip \ H \ 'ip) 
with the corresponding product states ( |3.l| ). For the states which minimize {ip \ 
H \ ip) we obtain the following results. If G < AE/2 there is only one mean-field 
ground state 

l^o) = |l)...|l), (3.6) 

with energy 

AF 

^0 - (V-o \H\^o)= -N—- (3.7) 



If G > AE/2 there are two degenerate mean-field ground states 

\^^) = |A^>...|A^) 

|^«) = |A«)...|A^), (3.9) 

where 



\^^) = |A^>...|A^) (3.8) 



1^') ^ 4-f i'>W^f i^> <"»> 



with energy 



|A«) = a-|A^), (3.11) 



A^ N fAE "■ ^ 
= -N 



^o--N—-^i—-G] . (3.12) 



By defining the critical value 



Gcr = ^, (3.13) 



we distinguish the following two cases. For G G (0, Gcr) the ground state of the 
system is approximated by a condensate of delocalized symmetric single particle 
states corresponding to the ground state of an isolated molecule. For G G {Gcr, oo) 
we have two different condensates which approximate the ground state of the sys- 
tem. The corresponding single particle states transform one into the other under 



the action of the inversion operator a^, as shown by Eq. (3.11), and, for G ^ Gc 
they become localized 



hm |A^) = \L) (3.14) 



AE 
G 



-0 



lim |A^) = \R). (3.15) 



Ag 
G 



+0 
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The above results suggest the existence of a quantum phase transition at G 



Qcr = -Tr-27^- (3-16) 



Gcr- According to Eq. (2.13), this phase transition can be obtained by increasing 
the density q of the gas above the critical value 

AS 1 

In the next Section we will discuss the behavior of the inversion line obtained in 
the range < g < Qcr when the gas, still being in the region characterized by 
delocalization, approaches the phase transition. 

4 Absorption spectrum 

When the gas is exposed to an electro-magnetic radiation of frequency vo = 
a;o/27r, we add to the Hamiltonian (2.10|) the perturbation 



Hem{t) = ef{t)S (4.1) 

where e is a small parameter, 

/(t)=0(t)cos(cc>oi), (4.2) 

d{t) being the Heaviside function, and S* is a time independent operator. In the 
long-wavelength approximation, Ao 3> d, where d is the molecule size, we have 

N 

S ^^li® ...®(jI®...®1n. (4.3) 

i=l 

At zero temperature the absorption spectrum induced by the electro-magnetic 
perturbation consists in a series of lines of frequency Vn = {En — EQ)h~^ correspond- 
ing to transitions of the system from the ground state of H to the n-th excited level. 
The corresponding transition probability is proportional to the square of the matrix 
element 

Mo„ = (Vo|^|V'«). (4.4) 

Since we are not able to determine exactly the eigenvalue problem i/Vn = Enipn for 
the N particle Hamiltonian ( ^.10| ), we will find an approximate mean- field solution 
in the framework of the linear response theory. 
Let us define the function 

5(t)^(V^(i)|5|V(i)>, (4.5) 

where |^(i)) = exp[— i(iJ + He,n)t/fi\\t(jo). By expanding S{t) in powers of e, for 
the first order term Si (t) we find iQ 

5l(c^) = />)7^(u;), (4.6) 

where Si{lu) and /(cj) are the Fourier transforms of 5i(t) and f{t) and 

7^(w) = J2 l^^o. 



1 1 



huj - {En - Eq) +i6 huj+ {En - Eq) + iS 



(4.7) 



with (5^0+. The frequencies of the spectroscopic lines and the associated tran- 
sition probabilities are therefore given by the poles and the residues of TZ{uj). In 
order to determine TZ{uj), we evaluate iSi(u;) in a mean-field approximation and 



then use (4.6). 
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The mean- field approximation for the time-dependent state {ipit)) is obtained 
by extremizing the action functional 

{^{t) I -ihdt + H + H,„, I ,lj{t))dt (4.8) 

/o 

in the class of states of the form 

m)) = \X{t))...\Xit)). (4.9) 

The single-particle state |A(i)) is thus determined by the nonlinear Schrodinger 
equation 

'^^^ = ihim) + e/W^1 |A(t)), (4.10) 



where h{X) is given by (3.5). Let us assume that the gas is in the quantum region 



G e (0, Gc r), so that the mean-field ground state is 1-00) = |1) ■ • • |1)- Since ["(/"(O)) 



iV'o): Eq- ( 4.10 ) has to be solved with the initial condition 

|A(0)) = |1). (4.11) 

The calculation shows that TZ{lo) has a single pole corresponding to a line in the 

(4.12) 
(4.13) 



absorption spectrum 


with frequency 

^= h V 








1 
1- 


2G 

AE 


The residue of 7^(w) 


at this pole is 


N 








s/^- 


2G 
AE 





Equation (f4.12| ) shows a dependence of u on the coupling constant G. For 
G = we have P = AEh^^, namely the inversion line frequency expected from the 
quantum theory of an isolated molecule. In the region < G < Gcr-, v decreases by 
increasing G and vanishes for G = Gcr- In figure || we show the described behavior 
of the inversion line frequency as a function of the gas density g = G/Cfi^. 



The residue (4.13) multiplied by the photon energy hv provides, up to a con- 



stant, the intensity / of the inversion line 

N 
I (xhD—^^^^ NAE. (4.14) 

Note that the divergence at G = Gcr shown by the residue is cancelled by the factor 
hD which vanishes in the same limit. 

5 Optical activity and superselection rules 

Our model does not distinguish molecules of the type XY3 from the substituted 
derivatives XYWZ . In fact, we believe that this distinction should not matter for 
the disappearance of the inversion line. However, an important difference between 
the two cases is that for XY^ the localized states can be obtained one from the 
other either by a rotation or by space inversion, while for XYWZ they can be 
connected only by space inversion. This implies that XYWZ molecules, whenever 
become localized, i.e. at a density greater than the critical value, are chiral and 
must show stable optical activity. Stable optical activity is commonly observed for 
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1.5 



Figure 1 Inversion-line frequency ;/ as a function of the gas density q. 



some substituted derivatives of AsH^ for which the critical density is exceedingly 
small. For the theory of optical activity we refer, for example, to po| . 

To explain the stability of the optical activity, i.e. why the localized states are 
stable, although they are not eigenstates of the single-molecule Hamiltonian, an 
environment induced superselection rule has been invoked [^, |lO|, |ll|, |l^ . 

Following ||ri|, Q we say that in a system with infinitely many degrees of 
freedom there is a superselection rule between two different sectors of a Hilbert 
space if there are neither spontaneous transitions between their state vectors, nor 
any measurable quantity with matrix element different from zero between them. 
This implies that any coherent superposition of states among the two sectors is 
physically meaningless and that we can construct a classical observable (i.e. an 
observable with zero dispersion on all the states) which acquires the same value 
on all the vectors of each sector. Examples are the charge and the Boson-Fermion 
superselection rules. The first refers to the impossibility of coherent superpositions 
of states with different charge, the second to the impossibility of superposing states 
with integer and half-integer angular momentum. 

With respect to this definition, in the limit of an infinite number of molecules 
(thermodynamic limit) our model has a superselection rule similar to that found by 
Pfcifer analyzing a single molecule interacting with the radiation field P|. Above 
the critical density, the Hilbert space separates into two sectors generated by 



the ground state vectors given in mean field approximation by (3.9). These sec- 
tors, which we call TYl and 'Hr, cannot be connected by any operator involving 
a finite number of degrees of freedom (local operator). Let us define chirality 
X = hm^v-^oo -^ Si=i li fH) •■• fX" erf (8) ... (8) In- In the limit A^ ^ oo we have 
xV" — iV" for 4' in Tii and TLr, respectively. Optical activity is expected when 
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From the above definition the concept of superselection rule is applicable only 
to pure states. A superselection rule, however, puts some restrictions in the con- 
struction of ensembles (density matrices) because these cannot contain mixed terms 
corresponding to vectors belonging to different sectors of the Hilbert space. Indeed, 
at finite temperature, we find that in the mean field approximation the free energy 
of our model is minimized by two density matrices which have vanishing matrix 
elements between the two sectors TYl and Tij^ [071. Any other admissible density 
matrix will be a convex combination of these two. 

A different point of view consists in considering, instead of the whole gas, a 
subsystem composed by a single molecule or a finite number of them still interacting 
with the rest of the gas. In this case a superselection rule in the Hilbert space of the 
subsystem can emerge dynamically. By this we mean that there exist projection 
operators such that the off-diagonal contributions of the subsystem density matrix 
vanish when time t —^ oo |2l[] . This alternative point of view, usually referred to 
as environment induced superselection rule |22f| , can be presumably implemented 
in our model if we describe the interaction of a molecule with the rest of the gas 
by means of a Lindblad equation psl . 



It is a pleasure to dedicate this paper to Sergio Doplicher on the occasion of his 
60th birthday and we thank him for a very interesting discussion on superselection 
rules. 
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